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Abstract 

Recent academic work has developed a method to determine, in real time, 
if a given stock is exhibiting a price bubble. Currently there is speculation in 
the financial press concerning the existence of a price bubble in the aftermath 
of the recent IPO of Linkedln. We analyze stock price tick data from the 
short lifetime of this stock through May 24, 2011, and we find that Linkedln 
has a price bubble. 

1 Introduction and Summary 

It has been well documented in the financial press that a methodology is needed 
that can identify an asset price bubble in real time. William Dudley, the President 
of the New York Federal Reserve, in an interview with Planet Money [3] stated 
"...what I am proposing is that we try to identify bubbles in real time, try to develop 
tools to address those bubbles, try to use those tools when appropriate to limit the 
size of those bubbles and, therefore, try to limit the damage when those bubbles 
burst." 

It is also widely recognized that this is not an easy task. Indeed, in 2009 the 
Federal Reserve Chairman Ben Bernanke said in Congressional Testimony [I] "It is 
extraordinarily difficult in real time to know if an asset price is appropriate or not" . 
Without a quantitative procedure, experts often have different opinions about the 
existence of price bubbles. A famous example is the oil price bubble of 2007/2008. 
Nobel prize winning economist Paul Krugman wrote in the New York Times that 
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it was not a bubble, and two days later Ben Stein wrote in the same paper that it 
was. 

Although not yet widely known by the finance industry, the authors have recently 
developed a procedure based on a sophisticated mathematical model for detecting 
asset price bubbles in real time (see We have successfully back-tested our 

methodology showing the existence of price bubbles in various stocks during the 
dot-com era of 2000 to 2002. We also showed that some stocks in that period that 
might have suspected of being bubbles, were in fact not. But, we have not yet tested 
our method on stocks in real time. That is, we have not tested them until now. 

Inspired by a New York Times article [2] discussing whether or not in the aftermath 
of the Linkedln IPO the stock price had a bubble, we obtained stock price tick 
data from Bloomberg. And, we used our methodology to test whether Linkedln's 
stock price is exhibiting a bubble. We have found, definitively, that there is a price 
bubble! 

Our method consists of assuming a general (and generally accepted) evolution for 
the stock price, estimating its volatility using state of the art estimators, and then 
extrapolating the volatility function to see if a certain calculus integral based on 
the volatility is finite or infinite. If it is finite, the stock price has a bubble; if the 
integral in question is infinite then the stock is not undergoing bubble pricing; our 
test is indeterminate in a small set of cases where the volatility tends to infinity at 
a rate between both possibilities. In the case of Linkedln, the volatility function is 
well inside the bubble region. There is no doubt about its existence. 

Our results can be summarized in the following graph showing an extrapolation of 
an estimated volatility function for Linkedln's stock price. Put simply, the theory 
developed in and [I] tells us that if the graph of the volatility versus the stock 

price tends to infinity at a faster rate than does the graph of f(x) = x, then we have 
a bubble. Below we have a graph of the volatility coefficient of Linkedln together 
with its extrapolation, and the reader can clearly see that the graph indicates the 
stock has a price bubble. 

The blue part of the graph is the estimated function, and the red part is its extrap- 
olation using the technique of Reproducing Kernel Hilbert Spaces (RKHS). 

These Linkedln results illustrate the usefulness of our new methodology for detecting 
bubbles in real time. Our methodology provides a solution to the problems stated 
by both Chairman Bernanke and President Dudley, and it is our hope that they will 
prove useful to regulators, policy makers, and industry participants alike. 

Acknowledgement: The authors thank Peter Carr and Arun Verma for help in 
obtaining quickly the tick data for the stock Linkedln. 
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Figure 1: Estimation and Extrapolation of the Volatility Function 

The Empirical Test 

In this section we first review the methodology contained in [H] to determine whether 
or not Linkedln's stock price is experiencing a bubble, and then we apply this 
methodology to Linkedln minute by minute stock price tick data obtained from 
Bloomberg. We conclude that Linedln's stock is indeed experiencing a price bub- 
ble. 

1.1 The Methodology 

The methodology is based on studying the characteristics of Linkedln's stock price 
process. If Linkedln's stock price is experiencing a bubble, it can be shown [6] 
that the stock price's volatility will be unusually large. Our empirical methodology 
validates whether or not Linkedln's stock volatility is large enough. 

To perform this validation, we start by assuming that the stock price follows a 
stochastic differential equation of a form that makes it a diffusion in an incomplete 
market (see ^ below). The assumed stock price evolution is very general and fits 
most stock price processes reasonably well. For identifying a price bubble, the key 
characteristic of this evolution is that the stock price's volatility cr(St) is a function 
of the level of the stock price. 

Next, we use sophisticated methods to estimate the volatility coefficient a(x). Since 
the data is necessarily finite, we can estimate the values of a(x) only on the set of 
stock prices observed (which is a compact subset of M + ) We use two methods to 
extend a to all of R+. One, we use parametric methods combined with a comparison 
theorem. Two, we use an indexed family of Reproducing Kernel Hilbert Spaces 
(RKHS) combined with an optimization over the index set to obtain the best possible 
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extension given the data (this is a sort of bootstrap procedure). 

This "knowledge of cr(St)" then enables us to determine whether the stock price 
is a martingale or is a strict local martingale under any of an infinite collection 
of risk neutral measures. If it is a strict local martingale for all of the risk neutral 
measures (which corresponds to a certain calculus integral being finite by Feller's test 
for explosions), then we can conclude that the stock price is undergoing a bubble. 
Otherwise, there is not a stock price bubble. 



1.2 The Estimation 

We assume that Linkedln's stock price is a diffusion of the form 

dS t = a(S t )dW t + b(S t ,Y t )dt (1) 
dY t = s(Y t )dB t + g(Y t )dt (2) 

where W and B are independent Brownian motions. This model permits that under 
the physical probability measure, the stock price can have a drift that depends on 
additional randomness, making the market incomplete (see [B]). Nevertheless for this 
family of models, S satisfies the following equation for every neutral measure: 

dS t = a(S t )dW t 

Under this evolution, the stock price exhibits a bubble if and only if 

x 

—, dx < oo for all a > 0. (3) 
a 2 (x) 

We test to see if this integral is finite or not. 

To perform this test, we obtained minute by minute stock price tick data for the 
4 business days 5/19/2011 to 5/24/2011 from Bloomberg. There are exactly 1535 
price observations in this data set. The time series plot of Linkedln's stock price is 
contained in Figure [2j The prices used are the open prices of each minute but the 
results are not sensitive to using open, high or lowest minute prices instead. 

The maximum stock price attained by Linkedln during this period is $120.74 and 
the minimum price was $81.24. As evidenced in this diagram, Linkedln experienced 
a dramatic price rise in its early trading. This suggests an unusually large stock 
price volatility over this short time period and perhaps a price bubble. 

Our bubble testing methodology first requires us to estimate the volatility function 
a using local time based non-parametric estimators. We use two such estimators. 
We compare the estimation results obtained using both Zmirou's estimator (see 
Theorem 1 in [J]) and the estimator developed in jl] (see Theorem 3 in the same 
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Figure 2: Linkedln Stock Prices from 5/19/2011 to 5/24/2011. (The observation 
interval is one minute.) 

reference). The implementation of these estimators requires a grid step h n tending 
to zero, such that nh n — > oo and nh^ — > for the former estimator, and nh? n — > oo 
for the later one. We choose the step size h n = -\ so that all of these conditions are 

simultaneously satisfied. This implies a grid of 7 points. The statistics are displayed 
in figure [3j 
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JKP Estimator 
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Figure 3: Non-parametric Volatility Estimates. 

Since the neighborhoods of the grid points $118,915 and $125,764 are either not vis- 
ited or visited only once, we do not have reliable estimates at these points. Therefore, 
we restrict ourselves to the grid containing only the first five points. We note that the 
last point in the new grid $112,065 still has only been visited very few times. 

When using Zmirou's estimator, confidence intervals are provided. The confidence 
intervals are quite wide. Given these observations, we apply our methodology twice. 
In the first test, we use a 5 point grid. In the second test, we remove the fifth 
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point where the estimation is uncertain and we use a 4 point grid instead. The 
graph in figure [4] plots the estimated volatilities for the grid points together with 
the confidence intervals. 
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Figure 4: Non parametric Volatility Estimation Results. 



The next step in our procedure is to interpolate the shape of the volatility function 
between these grid points. We use the estimations from our non parametric estima- 
tor with the 5 point grid case. For the volatility time scale, we let the 4 day time 
interval correspond to one unit of time. This scaling does not affect the conclusions 
of this paper. When interpolating one can use any reasonable method. We use both 
cubic splines and reproducing kernel Hilbert spaces as suggested in [4], subsection 
5.2.3 item (ii). The interpolated functions are in figure [5j 

From these,we select the kernel function K\ T as defined in Lemma 10 in [4J, and we 
choose the parameter r = 6. 




Figure 5: Interpolated Volatility using Cubic Splines and the RKHS Theory 



The next step is to extrapolate the interpolated function a b using the RKHS theory 
to the left and right stock price tails. Define f(x) = -^-^ and define the Hilbert 
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space 



H n = H n ([0,oo[) ={fe C n ([0,oo[) | lim x k f {k \x) = for all < k < n - 1} 

where n is the assumed degree of smoothness of /. We also need to define an inner 
product. A smooth reproducing kernel q RP (x,y) can be constructed explicitly (see 
Proposition 2 in [1]) via the choice 

y n f {n \y)y n g {n \y) d y 

n\ n\ w(y) 

where w(y) = -^ is an asymptotic weighting function. We consider the family 
of RKHS H nm = (H n , (,) n ,m), in which case the explicit form of is provided 
in Proposition 2 in [3] in terms of the Beta and the Gauss's hypergeometric func- 
tions. 

For n G {1, 2} fixed, we construct our extrapolation a = a m as in jl], 5.2.3 item (iv), 
by choosing the asymptotic weighting function parameter m such that f m = -\- is in 
H n , m , &m exactly matches the points obtained from the non parametric estimation, 
and a m is as close (in norm 2) to <j b on the last third of the bounded interval where 
a b is defined. Because of the observed kink and the obvious change in the rate of 
increase of a b at the forth point, we choose n — 1 in our numerical procedure. The 
result is shown in figure [6] 
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Figure 6: RKHS Based Extrapolation of a b 
We obtain m = 9.42. 

From Proposition 3 in [I], the asymptotic behavior of a is given by 

M 

lim x m+1 f{x) = n 2 B{m + 1, n) a 

x—tao ' * 

i=l 
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where M = 5 is the number of observations available, B is the Beta function, and 
the coefficients (ci)i<i<M are obtained by solving the system 

M 

J2d<ln^i,Xk) = f{xk) for all 1 < k < M 

i=l 

where (a?i)i<i<M is the grid of the non parametric estimation, f(xk) = a 2( Xk ) and 
a(xk) is the value at the grid point x^ obtained from the non-parametric estimation 
procedure. This implies that a is asymptotically equivalent to a function propor- 
tional to x a with a = "4p, that is a = 5.21. This value appears very large, however, 
the proportionality constant is also large. The q's are automatically adjusted to ex- 
actly match the input points (xi, /(a?i))i<i<M- 

We plot below the functions with different asymptotic weighting parameters m ob- 
tained using the RKHS extrapolation method, without optimization. All the func- 
tions exactly match the non-parametrically estimated points. 




Figure 7: Extrapolated Volatility Functions using Different Reproducing Kernels 

The asymptotic weighting function's parameter m = 9.42 obtained by optimization 
appears in figure [7] to be the estimate most consistent (within all the functions, 
in any Hilbert Space of the form H l m , that exactly match the input data) with a 
"natural" extension of the behavior of o b to IR + . The power a = 5.21 implies then 
that Linkedln stock price is currently exhibiting a bubble. 

Since there is a large standard error for the volatility estimate at the end point 
$112,065, we remove this point from the grid and repeat our procedure. Also, the 
rate of increase of the function between the last two last points appears large, and 
we do not want the volatility's behavior to follow solely from this fact. Hence, we 
check to see if we can conclude there is a price bubble based only on the first 4 
reliable observation points. We plot in figure [8] the function o~ b (in blue) and its 
extrapolation to M + , a (in red). 
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Figure 8: RKHS based Extrapolation of a b 

Now M — 4. With this new grid, we can assume a higher regularity n = 2 and we 
obtain, after optimization, m = 7.8543. This leads to the power a = 4.42715 for 
the asymptotic behavior of the volatility. Again, although this power appears to be 
high given the numerical values (xk, f(xk))i<k<^, the coefficients (Q)i<j<4 and hence 
the constant of proportionality are adjusted to exactly match the input points. The 
extrapolated function obtained is the most consistent (within all the functions, in 
any #2,mj that exactly match the input data) in terms of extending 'naturally' the 
behavior of o~ b to 1R + . Again, we can conclude that there is a stock price bubble. 
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